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Using first-principles calculations, we identify magic wavelengths A for the 2s — 2p and 2s — 3p tran- 
sitions in lithium. The ns and np atomic levels have the same ac Stark shifts at the corresponding 
magic wavelength, which facilitates state-insensitive optical cooling and trapping. Tune-out wave- 
lengths for which the ground-state frequency-dependent polarizability vanishes are also calculated. 
Differences of these wavelengths between 6 Li and 7 Li are reported. Our approach uses high-precision, 
relativistic all-order methods in which all single, double, and partial triple excitations of the Dirac- 
Fock wave functions are included to all orders of perturbation theory. Recommended values are 
provided for a large number of Li electric-dipole matrix elements. Static polarizabilities for the 2s, 
2p, 3s, 3p, and 3d levels are compared with other theory and experiment where available. Uncertain- 
ties of all recommended values are estimated. The magic wavelengths for the uv 2s — 3p transition 
are of particular interest for the production of a quantum gas of lithium [Duarte et al., Phys. Rev. 
A 84, 061406R (2011)]. 



I. INTRODUCTION 

The alkali atoms, with one electron outside a closed 
shell core, have an oscillator strength sum, £/ = 1, to be 
distributed among all optical transitions from the ground 
state at energies below the threshold for core excitation 
Most of this oscillator strength is concentrated in 
the lowest, ns — np, resonance transitions, where n = 2 
for Li, 3 for Na, etc. The higher resonance transitions, 
ns — n'p with n' > n, are thus usually much weaker than 
their conterparts in non-alkali atoms, and their natural 
linewidths are narrower. In some cases this is accompa- 
nied by anomalous intensity ratios of the fine-structure 
doublet lines, a phenomenon first explained by Fermi 

Recent experiments in 6 Li and 40 K Q have shown 
the advantages of laser cooling with higher resonance 
lines for reducing the temperature and increasing the 
phase space density of optically trapped alkali atoms. 
The narrow linewidths of these transitions reduce the 
Dopplcr cooling limit compared to that of the lowest res- 
onance line. Two direct tests were made using similar 
schemes: laser cooling of the gas in a magneto-optical 
trap (MOT) using the usual ns — > np 3 / 2 D2 transition 
(the wavelengths of which are 671 nm for 6 Li and 767 nm 
for 40 K ) , followed by transfer of the gas into a MOT op- 
crating at the ns — > (n + l)pz/2 UV line (respectively 
323 nm and 405 nm). In both cases, temperature reduc- 
tions by a factor of about five and phase-space density 
increases by at least a factor of ten were observed. 

In this paper we discuss aspects of this cooling scheme 
as applied to lithium. In order to be able to continue to 
laser cool on the UV transition when a dipole optical trap 
is turned off, the ac Stark shifts of the 2s and 3p 3 / 2 levels 
due to trap light have to be nearly the same, resulting in a 



sufficiently small differential ac Stark shift on the cooling 
transition to allow for efficient and uniform cooling 
The ac Stark shifts of these states are generally different 
and may be of different signs, leading to heating. The 
same problem, i.e. different Stark shifts of two states, af- 
fects optical frequency standards based on atoms trapped 
in optical lattices, because it can introduce a significant 
dependence of the measured frequency of the clock tran- 
sition upon the lattice wavelength. The idea of using 
a trapping laser tuned to a "magic" wavelength, A mag i c , 
at which the ac Stark shift of the clock transition van- 
ishes, was first proposed in Refs. [j| @. The use of the 
magic wavelengths is also important in order to trap and 
control neutral atoms inside high-Q cavities in the strong 
coupling regime with minimum decoherence for quantum 
computation and communication schemes 0]. 

The goal of the present work is to provide the list of 
precise magic wavelengths for Li UV 2s — 3pj transitions 
in convenient wavelength regions. We also provide a list 
of available magic wavelength for the 2s — 2pj transitions, 
calculate dc and ac polarizabilities for several low-lying 
states, and provide recommended values for a number 
of relevant electric-dipole transitions. We present results 
for both 6 Li and 7 Li to illustrate the possibilities for dif- 
ferential light shifts between the two isotopes. 

The magic wavelengths for a specific transition are lo- 
cated by calculating the ac polarizabilities of the lower 
and upper states and finding their crossing points. Magic 
wavelengths for lowest np — ns transitions in alkali-metal 
atoms from Na to Cs have been previously calculated in 
using a relativistic linearized coupled-cluster method. 
The data in Ref. 8] provided a wide range of magic wave- 
lengths for alkali-metal atoms. A bichromatic scheme 
for state-insensitive optical trapping of Rb atom was ex- 
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plored in Ref. Q. In the case of Rb, the magic wave- 
lengths associated with monochromatic trapping were 
sparse and relatively inconvenient. The bichromatic ap- 
proach yielded a number of promising magic wavelength 
pairs. 

Wc have also carried out calculations of the dc polariz- 
abilities of the 2s, 2pj, 3s, 3pj, and 3dj states to compare 
with available experimental and high-precision theoreti- 
cal values. We refer the reader to the recent review 
for extensive comparison of various results for Li polariz- 
abilities. Here, we provide comparison of our results with 
the most recent calculations carried out using Hylleraas 
basis functions [IlTfl5| since this approach is expected to 
produce the most accurate recommended values. Com- 
parison with the recent coupled-cluster calculation of flU ] 
and configuration interaction calculations with core po- 
larization (CICP) [HI is also included. 

We also provide values for the tune-out wavelengths, 
A ZC ro, for which the ground-state frequency-dependent 
polarizability of Li atoms vanishes. At these wavelengths, 
an atom experiences no dipole force and thus is unaf- 
fected by the presence of an optical lattice. These wave- 
lengths were first discussed by LeBlanc and Thywissen 
17j, and have recently been calculated for alkali- metal 
atoms from Li to Cs in 18]. Our present work provides 
accurate predictions of the tune-out wavelengths A zero for 
both 6 Li and 7 Li and explores the differences between the 
first tune-out wavelengths for these isotopes in detail. 

We start with the discussions of the calculation of 
electric-dipole matric elements, static and dynamic po- 
larizabilities as well as their uncertainties. The static 
values are compared with other theory and experiment 
in Section Mil The magic wavelengths are discussed in 
Section [TV] 



II. METHOD 

The background to our approach to calculation of 
atomic polarizabilities was discussed in Refs. j8l [l8l. [l9| . 
Here, we summarize points salient to the present work. 
The frequency-dependent scalar polarizability, a(w), of 
an alkali-metal atom in its ground state v may be sepa- 
rated into a contribution from the core electrons, a corc , a 
core modification due to the valence electron, a vc , and a 
contribution from the valence electron, a v (ui). The core 
polarizability depends weakly on to for the frequencies 
treated here, since core electrons have excitation energies 
in the far-ultraviolet region of the spectrum. Therefore, 
we approximate the core polarizability by its dc value as 
calculated in the random-phase approximation (RPA). 
The accuracy of the RPA approach has been discussed in 
[Tol | ■ T ne core polarizability is corrected for Pauli block- 
ing of core-valence excitations by introducing an extra 
term a vc . For consistency, this is also calculated in RPA. 
The valence contribution to frequency-dependent scalar 
ao and tensor oti polarizabilities is evaluated as the sum 
over intermediate k states allowed by the electric-dipole 
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and (k \\d\\ v) are the reduced electric-dipole matrix ele- 
ments. In these equations, oj is assumed to be at least 
several linewidths off resonance with the corresponding 
transitions. Linear polarization is assumed in all calcu- 
lation. 

Unless stated otherwise, we use the conventional sys- 



tem of atomic units, a.u., in which 
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reduced Planck constant fi have the numerical value 1. 
Polarizability in a.u. has the dimension of volume, and 
its numerical values presented here are expressed in units 
of (Iq, where ao ~ 0.052918 nm is the Bohr radius. 
The atomic units for a can be converted to SI units via 
a/h [Hz/(V/m) 2 ]=2.48832xl0- 8 a [a.u.], where the con- 
version coefficient is Aite^a^/h and the Planck constant h 
is factored out. 

We use the linearized version of the coupled clus- 
ter approach (also refereed to as the all-order method), 
which sums infinite sets of many-body perturbation the- 
ory terms, for k with principal quantum number n < 26. 
The 2s — np, 2p — nl, 3s — nl, Zp — nl, and 3d — nl tran- 
sitions with n < 26 are calculated using this approach. 
Detailed description of the all-order method is given in 
Refs. |3, HH. 

We use experimental values of energy levels up to 
n = 12 taken from [13, HH, H3] and theoretical all-order 
energy levels for n — 13 — 26. The remaining contribu- 
tions with n > 26 are calculated in the Dirac-Fock (DF) 
approximation. We use a complete set of DF wave func- 
tions on a nonlinear grid generated using B-splines con- 
strained to a spherical cavity. A cavity radius of 220 ao is 
chosen to accommodate all valence orbitals with n < 13 
so we can use experimental energies for these states. The 
basis set consists of 70 splines of order 11 for each value 
of the relativistic angular quantum number k. 

The evaluation of the uncertainty of the matrix ele- 
ments in this approach was described in detail in [2l| . 
Four all-order calculations were carried out. Two of these 
were ab initio all-order calculations with and without the 
inclusion of the partial triple excitations. Two other cal- 
culations included semiempirical estimate of high-order 
correlation corrections starting from both ab initio runs. 
The spread of these four values for each transition de- 
fines the estimated uncertainty in the final results. Since 
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TABLE I: Absolute values of the reduced electric-dipole ma- 
trix elements in Li and their uncertainties. Units: a.u. 
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high-order corrections are small for some Li transitions, 
in particular for the excited states, the uncertainty esti- 
mate for the lower transition was used for the other tran- 
sitions. For example, the 2p — 3s uncertainty was used for 
the other 2p — ns transitions with n > 3. This procedure 
insures that all uncertainty estimates are larger than the 
expected numerical accuracy of the calculations. The un- 
certainties for the small 2s — npj and 3s — n\Pj matrix 



elements with n > 2 and n\ > 3 were estimated as 10% of 
the total correlation correction. For these transitions, the 
procedure used to estimate high-order corrections does 
not estimate all dominant contributions. Placing the un- 
certainty estimate at 10% of the correlation correction 
for these transitions ensures that all missing correlation 
effects do not exceed this rough uncertainty estimate. 
The matrix element calculations have been carried out 
for a fixed nucleus using a Fermi distribution for a nuclear 
charge. These matrix elements are used in all polarizabil- 
ity calculations in the present work. We have estimated 
that the 6 Li - 7 Li isotope shift correction to the matrix 
elements is well below our estimated uncertainties. 

The absolute values of the reduced electric-dipole Li 
matrix elements used in our subsequent calculations and 
their uncertainties are listed in a.u. in Table fl] We note 
that we list only the most important subset of the ma- 
trix elements. We have calculated a total of 474 matrix 
elements for this work. 



TABLE II: Comparison of the present values of static scalar 
ao and tensor 02 7 Li polarizabilities with other calculations 
and experiment. The final entry is the difference of the 2p!/ 2 
and 2s polarizabilities for which there is a precise experimen- 
tal determination using dc Stark shift measurement. Units: 
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III. LI POLARIZABILITIES 

We start with a brief overview of the recent calcula- 
tions of Li polarizabilities. The highest accuracy attained 
in ab initio atomic structure calculations is achieved by 
the use of basis functions which explicitly include the in- 
terelectronic coordinate. Difficulties with performing the 
accompanying multi-center integrals have effectively pre- 
cluded the use of such basis functions for systems with 
more than three electrons. Correlated basis calculations 
are possible for lithium since it only has three electrons. 
Consequently it has been possible to calculate Li polar- 
izabilities to very high precision [TTI4T31 ITsT ] . The most 
accurate calculation of polarizability of lithium ground 
state ao(2s) = 164.1125(5) a.u. was carried out in [l3| . 
This calculation included relativistic and quantum elec- 
trodynamics corrections. The small uncertainty comes 
from the approximate treatment of quantum electrody- 
namics corrections. This theoretical result can be consid- 
ered as a benchmark for more general atomic structure 
methods and may serve as a reference value for the rel- 
ative measurement of polarizabilities of the other alkali- 
metal atoms [l3|, H3] ■ The uncertainty in the experimen- 
tal value of the polarizability 164.2(11) a.u. [24J is sub- 
stantially higher. The most stringent experimental test of 
Li polarizability calculations is presently the Stark shift 
measurement of the 2s-2p 1 / 2 transition by Hunter et al. 
25], which gave a polarizability difference of —37.15(2) 
a.u. Our coupled-cluster result —37.19(7) a.u. is in excel- 
lent agreement with the experimental polarizability dif- 
ference. 

The electric dipole polarizabilities and hyperpolariz- 
abilities for the lithium isotopes 6 Li and 7 Li in the 2s, 
2p, and 3d excited states were calculated nonrelativisti- 
cally using the variational method with a Hylleraas basis 
[TT| . The calculated 3d polarizabilities were found to be 
in significant disagreement with 2003 measurements [26[ . 
In order resolve this discrepancy, we have calculated 3d, 
static scalar and tensor polarizabilities. Our results are 
in excellent agreement with the Hylleraas calculations of 
11] indicating a problem with measurements reported in 
261 ] . The dc and ac dipole polarizabilities for Li atoms 
in the 2 s and 2p states were calculated using the same 
variational method with a Hylleraas basis in 12]. Cor- 
rections due to relativistic effects were also estimated to 
provide recommended values. 

Since the Hylleraas basis set calculations can not be 
carried out for larger system at the present time, Li pro- 
vides an excellent benchmark test of our coupled-cluster 
approach, and our procedures to estimate the uncer- 
tainties. Our approach is intrinsically relativistic unlike 
nonrelativistic Hylleraas calculations that requires sub- 
sequent estimate of the relativistic corrections. We note 
that the coupled-cluster method used in this work is not 
expected to provide accuracy that is competitive with 
what is ultimately possible to achieve with Hylleraas ba- 
sis for dc polarizability results. The accuracy of the 
coupled-cluster method was recently discussed in [28| . 
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FIG. 1: (Color online) The frequency-dependent polarizabili- 
ties of the Li 2s and 2p 1 / 2 states. The magic wavelengths are 
marked with arrows. 
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FIG. 2: (Color online) The frequency-dependent polarizabili- 
ties of the Li 2s and 2p 3 / 2 states. The magic wavelengths are 
marked with arrows. 



However, the coupled-cluster approach is applicable to 
more complicated systems, so comparison with Hylleraas 
calculations in Li provides stringiest tests of the method- 
ology as well as of numerical accuracy. Moreover, com- 
parison with Hylleraas benchmarks provides additional 
confidence in our estimates of the uncertainties in the 
values of magic wavelengths. Comparison of our results 
with most recent correlated basis set calculations 
is given in Table |TTJ All data are given for 7 Li since the 
differences between 6 Li and 7 Li values are smaller than 
our uncertainties. Our results are in excellent agreement 
with all Hylleraas values. 

Li electric-dipole and quadrupole polarizabilities, and 
van der Waals coefficients were calculated using relativis- 
tic coupled-cluster method with single, double, and par- 
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FIG. 4: (Color online) The frequency-dependent polarizabili- 
ties of the Li 2s and 3p 3 / 2 states. The magic wavelengths are 
marked with arrows. 



tial triple excitations, CCSD(T), in Ref. [29]. However, 
a mismatch of phases between the numerical and ana- 
lytical orbitals caused significant errors in the reported 
values. The later revision of this work [16j . where this 
problem was corrected, still yielded results that are in 
substantially poorer agreement with Hylleraas data than 
the values in the present work. The importance of using 
very accurate basis sets for coupled-cluster calculations 
of polarizabilities was discussed in a number of publica- 
tions including review [Kj . In our work, very large basis 
set with 70 orbitals for each partial wave was used result- 
ing in a very high numerical accuracy. Table |TT] includes 
the comparison with the revised CCSD(T) calculation of 
[l6j . and configuration interaction calculations with core 
polarization (CICP) [Hj]. 



TABLE III: Magic wavelengths for the 2s — npj transitions in 
6 Li and 7 Li. The 400 - 950 nm and 950 - 1300 wavelength 
ranges were considered for the 2s—2pj and 2s— 3pj transitions, 
respectively. The corresponding polarizabilities are given in 
a.u. 
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1107.793(5) 


1107.783(5) 


258.12(7) 






1287.25(6) 


1287.24(6) 


224.66(7) 


2s - 


3P3/2 


977.265(7) 


977.264(7) 


308.75(9) 


\ m j\ 


= 1/2 


998.826(3) 


998.823(3) 


297.50(8) 






1070.73(2) 


1070.72(2) 


269.01(8) 






1112.111(9) 


1112.102(8) 


256.97(7) 






1288.15(4) 


1288.15(4) 


224.55(7) 


2s - 


3P3/2 


976.240(7) 


976.239(7) 


309.33(9) 


\ m j\ 


= 3/2 


1068.76(2) 


1068.74(2) 


269.65(8) 






1285.88(6) 


1285.87(6) 


224.84(7) 



IV. MAGIC WAVELENGTHS 

We define the magic wavelength A mag i c as the wave- 
length for which the ac polarizabilities of two states in- 
volved in the atomic transition are the same, leading 
to vanishing ac Stark shift of that transition. For the 
ns — np transitions, a magic wavelength is represented 
by the point at which two curves, a ns (X) and a np (A), 
intersect as a function of the wavelength A. 

The total polarizability of the np 3 / 2 state depends 
upon its rrij quantum number. Therefore, the magic 
wavelengths need to be determined separately for the 
cases with rrij — ±1/2 and rrij — ±3/2 for the ns — np 3 / 2 
transitions, owing to the presence of the tensor con- 
tribution to the total polarizability of the np 3 / 2 state. 
The total polarizability for the np 3 / 2 states is given by 
a = ao — ot2 for rrij = ±1/2 and a = ao + «2 for the 
rrij = ±3/2 case. The uncertainties in the values of magic 
wavelengths are found as the maximum differences be- 
tween the central value and the crossings of the a ns ±5a ns 
and a np ±5a np curves, where the 5a are the uncertainties 
in the corresponding ns and np polarizability values. All 
calculation are carried out for linear polarization. The 
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TABLE IV: Tune-out wavelengths A zcro for 6 Li and 7 Li. The 



resonant wavelengths A rcs 


for relevant 


transitions are also 


listed. 


The wavelengths (in 


vacuum) are 


given in nm. 


Atom 


Resonance 




■^zcro 


6 Li 


2s — 2di /o 


670.992478 








670 987445(1) 




2s - 2p 3/2 


670.977380 










324.19(2) 




2s - 3pi/ 2 


323.3622 






2s - 3p 3/2 


323.361168 








274.920(10) 




2s - 4p 1/2 


274.2035 




7 Li 


2,s — 2?9i /o 


670.976658 








670 971625f2) 




2s - 2p 3/2 


670.961561 








324.19(2) 




2s - 3f> 1/2 


323.3572 






2s - 3p 3/2 


323.3562 








274.916(10) 




2s - 4p 1/2 


274.1998 





frequency-dependent polarizabilities of the Li ground and 
2pi/2 states for A = 500 — 900 nm are plotted in Fig. [T] 
The positions of the resonances are indicated by verti- 
cal lines with small arrows on top of the graph. There 
are only 3 resonances (2p — 3d, 2s — 2p, and 2p — 3s) in 
this wavelength region, resulting in only two magic wave- 
lengths above 500 nm for the 2s — 2pi/ 2 transition. These 
occur at 549.41(6) nm and 872.56(9) nm. The magic 
wavelengths are marked with arrows on all graphs. The 
respective values of the polarizabilities are also given for 
all magic wavelengths illustrated in this work. There are 
no other magic wavelengths for A > 900 nm since there 
are no more resonance contributions to both 2 s and 2p 1 / 2 
polarizabilities in this region. Both polarizabilities slowly 
decrease for A > 900 nm to approach their static value 
in the u) — > limit. 

The frequency-dependent polarizabilities of the Li 
ground and 2p 3 / 2 states for A = 500 — 900 nm are plotted 
in Fig. [21 The same designations are used in all figures. 
The magic wavelengths for the 2s — 2p 3 / 2 \mj\ = 1/2 
transition are similar to the ones for the 2s — 2p x / 2 tran- 
sition. In the \rrij\ = 3/2 case, there is no magic wave- 
length above 550 nm since there is no contribution from 
the 2s — 2p resonance. 

The magic wavelengths for the 2s — ip\/ 2 and 2s — ip^/2 
transitions for A = 950 — 1200 nm are illustrated in Fig. [3] 
and Fig. [4] All resonances indicated on top of the figures 
refer to the 3p — nl transitions, i.e. arrow labelled 7s in- 
dicates the position of the 3p — 7s resonance. A number 
of the magic wavelengths are available for the 2s — 3p 
uv transitions owing to a large number of resonance con- 
tributions to the 2>p polarizabilities in this region. The 
3s— nl resonances with n > 7 will yield other magic wave- 
lengths to the left of the plotted region, while the 3s — nl 



resonances with n < 6 will yield other magic wavelengths 
to the right of the plotted region. 

We have also calculated the ac polarizability of the 3s 
state for A = 950 — 1200 nm. We find that 3s polariz- 
ability is negative and large («-1000 a.u.) in this entire 
wavelength region with the exception of the very narrow 
wavelength interval between 1079.48 - 1079.50 nm (due 
to 3s — Apj resonances near 1079.49 nm). 

Table |TlT] presents the magic wavelengths for the 2s — 
2pj transitions in the 400 — 950 nm region, and the 
magic wavelengths for the 2s — 3pj transitions in the 
950 — 1300 nm region. We have carried out separate 
calculations for 6 Li and 7 Li by using the experimental 
energies for each isotope from [22, Ho] ■ The same theo- 
retical matrix elements were used for both isotopes, since 
the isotopic dependence of the matrix elements is much 
less than out uncertainty. We find that the differences 
between 6 Li and 7 Li magic wavelengths is very small and 
is smaller than our estimated uncertainties for almost 
all of the cases. We list both sets of data for the illus- 
tration of the IS differences. The polarizabilities at the 
magic wavelengths for 6 Li and Li are the same within 
the listed uncertainties, so only one set of polarizability 
data is listed. 

While our calculations are not sensitive enough to sig- 
nificantly differentiate between magic wavelengths for 6 Li 
and 7 Li, our values of the first tune-out wavelength are 
significantly different for the two isotopes [3]. There- 
fore, we consider Li tune-out wavelengths in more detail 
in the next section. 



V. LI TUNE-OUT WAVELENGTHS 

A tune-out wavelength A zcr o for a given state is one 
for which the ac polarizability of that state vanishes. In 
practice, we calculate ao(w) for a range of frequencies in 
the vicinity of relevant resonances and look for changes 
in sign of the polarizability within a given range. In the 
vicinity of the sign change, we vary the frequency until 
the residual polarizability is smaller than our uncertainty. 

In previous work [18] , we presented the tune-out wave- 
lengths for the ground state of Li. Here, we have calcu- 
lated the matrix elements to somewhat higher precision 
and have carried out more extensive study of their un- 
certainties. Therefore, we reevaluate the tune-out wave- 
lengths for Li in the present work. We have carried out 
separate calculations for 6 Li and 7 Li by using the exper- 
imental energies for each isotope from [H t3?| but the 
same matrix elements. Only the values of the first tune- 
out wavelengths are significantly affected by the IS. We 
illustrate this case separately in Fig. [5j The isotope shift 
of the 2s — 2p line is slightly larger than its fine structure. 
Therefore, the tune-out wavelength for one isotope cor- 
responds to a very deep trap for the other. We note that 
the first tune-out wavelength is very close to the Dl, D2 
resonances which would be a source of strong light scat- 
tering in an experiment. Applications of the tune-out 
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FIG. 5: (Color online) The frequency-dependent polarizabil- 
ities of the ground state of 6 Li and 7 Li. The tune-out wave- 
lengths are marked with arrows. 



wavelengths to sympathetic cooling and precision mea- 
surement were discussed in 1181. 



VI. CONCLUSION 

We have calculated the ground ns state and np state 
ac polarizabilities in Li using the relativistic linearized 
coupled-cluster method and evaluated the uncertainties 
of these values. The static polarizability values were 
found to be in excellent agreement with experiment and 
high-precision theoretical calculations with correlated ba- 
sis functions. We have used our calculations to identify 
the magic wavelengths for the 2s — 2p and 2s — 3p tran- 
sitions relevant to the use of ultraviolet resonance lines 
for laser cooling of ultracold gases with high phase-space 
densities. 
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